Great Basin Naturalist 56(3), © 1996, pp. 272-275 


INDICATORS OF RED SQUIRREL (TAMIASCIURUS HUDSONICUS) 
ABUNDANCE IN THE WHITEBARK PINE ZONE 


David J. Mattson! and Daniel P. Reinhart? 


ABSTRACT.—We investigated occupied squirrel middens and squirrel sightings and vocalizations as indicators of red 


squirrel (Tamiasciurus hudsonicus) abundance in the high-elevation whitebark pine (Pinus albicaulis) zone. Data were col- 


lected 1984-1959 from line transects located on 2 study sites in the Yellowstone ecosystem. We evaluated the performance 


of each measure on the basis of precision and biological considerations. We concluded that, of the 3 measures, active mid- 


dens were the best indicator of red squirrel abundance. We also observed that the density of active middens dropped by 


48%-66% between 1987 and 1989, following a severe drought and extensive wildfires that bummed one of the study sites 


during 19SS. 
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Whitebark pine (Pinus albicaulis) seeds are 
an important bear food that affects the survival 
and feeundity of grizzly bears (Ursus arctos) in 
the Yellowstone ecosystem. Use of pine seeds 
by grizzlies is almost entirely contingent upon 
the availability of cones cached in middens (i.e., 
larder hoards) by red squirrels (Tamiasciurus 
hudsonicus). Management of whitebark pine 
habitats for grizzlies has thus become contin- 
gent upon management of red squirrel popula- 
tions (Mattson and Reinhart 1994). 

We studied red squirrels in the whitebark 
pine zone using data collected from line tran- 
sects. Because these data included counts of 
middens, animals, and voealizations, we were 
able to evaluate the relative efficacy of these 3 
indicators of squirrel presence. We were inter- 
ested in identifving a “well-behaved” and rele- 
vant indicator of density to facilitate our inves- 
tigation of relationships between squirrel 
abundance and environmental factors such as 
midden use by grizzly bears. We were also 
interested in providing managers with an 
approach they could use to indicate squirrel 
abundance, short of using intensive methods 
that relied upon marked animals. 


STUDY AREA 
Our study area consisted of 2 sites, one 


loeated on the Mt. Washburn massif in north 
central Yellowstone National Park (44°47'N), 


and the other near Cooke City, Montana, 
immediately northeast of the park (45°00'N). 
These sites spanned the whitebark pine zone, 
from 2360 to 2870 m elevation. The whitebark 
pine zone borders upper timberline and is 
accordingly cold (average annual temperatures 
<0°C), often windy, and subject to deep (1-2 m) 
winter snow accumulations (Weaver 1990). 


MATERIALS AND METHODS 


Broad study objectives affected our transect 
design. We mapped the study area by habitat 
type-cover type strata based upon ground- 
truthed interpretation of 1:20,000 aerial pho- 
tography. The result was a fine-scale mosaic, 
with individual map polygons (forest stands) 
sometimes as small as individual squirrel terri- 
tories. To minimize effects of edge between 
different habitat types, we placed transects so 
as to maximize the number of right-angle inter- 
sections with stand boundaries as well as the 
amount of intersection with stand interiors. 
Because of this consideration and because for- 
est and meadow were variously intermixed, 
transect lines were of unequal length. 

We surveyed transects in the same order 
each vear, beginning after 10 August and end- 
ing prior to 28 September. Two observers walked 
permanently marked transect lines, with one ob- 
server primarily responsible for observations and 
the other primarily responsible for recording 
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data and keeping on line. At least one of the 
observers (the junior author) was the same 
during all years at both sites. Vocalizations and 
observed squirrels or middens (both active and 
inactive, by criteria of Finley [1969]) were 
recorded along with their estimated perpen- 
dicular distances from the transect line. Indi- 
vidual cone caches were not considered to be a 
“midden’ and were easily distinguished from 
these larger, more permanent features. 

We used the computer program TRANSECT 
(Burnham et al. 1980) to estimate densities. 
Individual transects constituted sample units 
for density calculations. As recommended by 
Burnham et al. (1980), we used the Fourier 
series, with 1—4 terms, to estimate distanee-to- 
line probability detection functions (g[x]). The 
distance at which we specified the limits of 
detectability for our measures (i.e., the eut-point) 
exerted considerable influence on the fit of the 
Fourier funetion to the observed detection dis- 
tribution. Accordingly, we varied cut-points to 
achieve the best fit to each vear- or site-spe- 
cific data set. Because data were collected 
from only 35 transects on the Mt. Washburn 
massif during 1984 and from 15 transects in 
the Cooke City area during 1984 and 1989 (com- 
pared to 57 and 21 transects, respectively, for 
all other years), we also calculated densities 
solely from these original 35 and 15 transects 
for all years so as to allow comparison with 
results from 1984. 


RESULTS 


We sampled the study area 5 yr, 1984-1987 
and 1989. During 1988 wildfires burned 562,000 
ha of the Yellowstone area, including 52% of 
the Mt. Washburn transects (47% severely). 
Transects on the Mt. Washburn area totaled 
18.9 km during 1984 (mean transect length [X] 
= 539 + 245 m [s]) and 29.8 km during the 
remaining 4 yr (X = 523 + 258 m). Similarly, 
during 1984 and 1989 transects on the Cooke 
City area totaled 16.4 kin (X = 1091 + 427 m) 
and 21.1 km during the remaining years (X = 
1005 + 405 m). We recorded 124 squirrel sight- 
ings, 641 voealizations, and 300 active middens 
on the Mt. Washburn study site and 54 sight- 
ings, 528 vocalizations, and 201 active middens 
on the Cooke City study site during the 5 study 
vears. The small number of sightings from the 
Cooke City site prevented us from estimating 
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annual densities from this measure for this 
aLOz, 

Total distance-to-line frequency distribu- 
tions for each of the 3 measures did not differ 
between the Mt. Washburn and Cooke City 
study sites (Mantel-Haenszel y2 for ordinal 
categories, P = (0.51, 0.61, and 0.35 for active 
middens, voealizations, and sightings, respec- 
tively). Perpendicular distributions of sightings 
and active middens peaked in the nearest (<10 
in) distance category, although the distribution 
of sightings more closely resembled a negative 
exponential and the distribution of middens a 
negative sigmoidal function. The majority (65% 
and 78%, respectively, by year and study site) 
of both these distributions were adcquately fit 
(y? test, P > 0.10) by a single-term Fourier 
funetion. Distributions of vocalizations peaked 
in the 2nd (11-20 m) distance category and 
were chuaracteristically (94%) fit by a 2- or 
higher-term Fourier fiunetion. In 3 (18%) in- 
stances we could not achieve an adequate fit 
by any model. 

Relationships among annual density esti- 
mates from the 3 measures were varied (Fig. 
1). On the Mt. Washburn site, mean sighting 
and vocalization densities were weakly corre- 
lated (r = 0.722), but tended to have overlap- 
ping 95% confidence intervals. Only 2 of 9 con- 
fidence intervals for the observed estimates (all 
years, for both the 1984 and inclusive samples) 
did not contain the line describing perfect cor- 
respondence (Fig. Id). In all but a single in- 
stance (Cooke City, 1984), mean midden densi- 
ties were greater than mean densities of the 
other 2 measures and were more strongly cor- 
related with sightings than vocalizations (r = 
0.981 versus r = 0.831, respectively, for tran- 
sects 1-57, Mt. Washburn; Fig. Ic). However, 
in this case, only two of nine 95% confidence 
intervals for midden and sighting densities in- 
cluded the possibility of perfect correspondence. 


CONCLUSIONS 


From these results we concluded that densi- 
ties calculated from active middens were more 
useful than densities calculated from the other 
2 measures for indicating red squirrel abun- 
dance. Our conclusion followed from the greater 
apparent detectability of middens compared to 
the squirrels themselves, the consistency with 
which a single-term Fourier function described 
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Fig. L. Relationships between annual estimates of density for active middens compared to vocalizations, (a) for Mt. 
Washburn and (b) for Cooke City, (e) sightings compared to active middens for Mt. Washburn, and (d) sightings compared 
to vocalizations, Mt. Washburn, 1984-1987 and 1989. Error bars correspond to 95% confidence intervals, solid circles to 
results from all transects, and open circles to results from the fewer transects established and first surveyed in 195-4. Diag- 


onals represent perfect correspondence between estimates. 


the probability detection distribution for mid- 
dens, and the resulting consistently smaller 
standard errors for the density estimates. In 
addition, scatter plots showed that aetive mid- 
den densities tended to be >0 when sighting 
and voealization densities were not. By impli- 
cation, vocalization and espeeially sighting 
densities were more likely to underestimate 
true squirrel densities; ie., at the same time 
that active middens clearly indicated the pres- 


ence of squirrels, sightings and voealizations 
could suggest there were none. 

Beeause red squirrel middens are nonmobile, 
often numerous, relatively easily observed, and 
typically associated with only one squirrel (Kil- 
ham 1954, M. Smith 1968, Wolff and Zasada 
1975, Vahle and Patton 1983), they are logical 
indicators of squirrel abundance. Furthermore, 
they do not suffer from sampling problems 
associated with weather, season, and time of day 
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TABLE 1. Estimated mean (n han!) and standard error (s¢) for densities of active middens on the VIL Washburn and 
Cooke City study sites, 1984-1987 and 1989, pereent coe Acie ait of variations for annual variation LOS4—L987. and percent 
decline in density from 1987 to 1989. Results are given for the transects established and surveyed during 19S. (1-35 and 
1-15) and for the larger sample of transeets surveyed during all other years, except for 1989, in the Cooke City area (1-57 


mo 1—21). 











Mt. Washburn Cooke City 
Trans. 1-35 Trans. 1-57 Trans. 1-15 ‘Trice | =21 

Yeur Mean (sz) \lean (sz) Mean (Se) \lean ise) 
19S4 Q.447 (0.083) — 0.498 (0.077) — 
1985 50% (0.110) 0.632 (0.08-1) 0.682 (0.116) 0.635 ((). O95) 
1986 ().219 (0.042) 0.426 (0.075) 0.5-18 (0.126) 0.5-10 ().095) 
1987 0.453 (0.093) 0.835 (0.143) 0.5 Fl (0.103) 0.790 (O.170) 
1989 ().234 (0.062) 0.285 (0.09%) 0.262 (Q. 137) — 
OV 

19S8--1987 30.4 32.6 18.9 19.2 
% decline 

1987-1989 48.3 66.0 Deo — 
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